Real Time Radiative Corrections to Charged Particle Decay Laws 



A. Widom and J. Swain 
Physics Department, Northeastern University, Boston MA U.S.A. 

Y.N. Srivastava 
Physics Department & INFN, University of Perugia, Perugia, Italy 



The real time exponential decay laws for meta-stable 
charged particles are shown to require radiative corrections. 
The methods employed are well known to be valid for radia- 
tively correcting Breit-Wigner line shapes. Radiative correc- 
tions contribute substantially to precision life time measure- 
ments of muons and pions when initially stopped in condensed 
matter. 
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I. INTRODUCTION 

It is well known p] that the exponential decay law Q 
for the survival probability Po(t) of a meta-stable parti- 
cle, 



P (t) « exp(-R), 



(1) 



arises from an energy distribution |3| which is nearly 
Lorentzian in shape 



dWo(E) 



dE 



r 



2W \(E-E ) 2 + (T/2y 



(2) 



i.e. 



P (t)= J e- tEt/n dW Q {E) 2 (exponential decay). (3) 

It is also very well known that when charged particles 
are involved in a decay process, the Lorentzian energy 
distribution dWo(E) of the meta-stable particle must be 
radiatively corrected (|] || to a new energy distri- 
bution dW(E). The new distribution has a much more 
"skewed" line shape than that of a simple Lorentzian. 
It then follows that the survival probability also has a 
radiative correction in real time; i.e. 



Pit) 



-iEt 



/n dW{E) 



(radiatively corrected) . 



(4) 



Our purpose is to discuss the real time consequences of 
radiation in the decay of meta-stable charged particles. 

Previous to this work, the notion that radiative cor- 
rections have important implications (7| for observations 
of quantum noise in a-decays and /3-decays in heavy nu- 
clei was theoretically developed M for real time counting 



rates. Strong experimental evidence for quantum noise 
in nuclear /3-decay counting rates has been reported j^] 

The time scales of quantum noise observations in nu- 
clear physics are of the order of a few hours to a few days. 
These long times are still much less than the nuclear life 
times, which in turn are of the order of a few years. For 
the application of radiative corrections to real time decay 
laws in high energy physics, e.g. for the weak decay 



/' 



or the weak decay 



/' 



(5) 



(G) 



the available absolute measurement time scales are con- 
siderably reduced. Thus, the methods of detecting real 
time radiative corrections should be somewhat different 
from those methods used for nuclear weak decays. 

In Sec. II, the radiatively corrected survival probability 
P(t) will be calculated from the mean number of soft 
photons dN(uj) radiated into the bandwidth dio. The 
properties of the resulting distribution are discussed in 
Sec. Ill, in terms of the time dependent transition rate 
for decay "f(t), defined as 



P{t) = exp ( - / j(s)ds 



(7) 



In Sec. Ill we illustrate the computation of the time de- 
pendent transition rate for a simple two body decay of 
a charged particle, e.g. Eq.(5). In Sec. IV we discuss 
the soft photon emission coupling strength as a func- 
tion of photon frequency and the velocity of the pro- 
duced charged particle both for insulators and conduc- 
tors. We discuss the soft photon emission in both the 
Bremsstrahlung and Cerenkov regimes. In Sec. V, we 
find (for times short on the scale of inverse soft photon 
frequencies) that there exists a "hot spot" in the decay 
rate; i.e. the decay rate per unit time j(t) has a sharp 
peak for short times. In the long time limit, the decay 
rate settles down to T, which determines the intrinsic life 
time. In the vacuum, as well as in materials, the de- 
cay rates exhibit a long time inverse t asymptotic time 
dependence which may be described by 



7(*) 



as t 



CO. 



(8) 
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In Eq.(8), the over all coupling strength Pq for the soft 
photons may be computed from the mean number of ra- 
diated photons dN(u) in a bandwidth du; It is 



f3 = hm w[ — 

lu^O \ WU 



(9) 



In the concluding Sec. VI, we discuss why the notion of 
radiative corrections to real time decay measurements in 
high energy physics is surely worthy of further experi- 
mental study. 



II. ENERGY DISTRIBUTIONS AND SURVIVAL 
PROBABILITIES 

Let VP denote the internal wave function of an unstable 
charged particle in the center of mass frame. The energy 
distribution of the state is given by 

dW(E) = (V,6(E- H)m)dE. (10) 

The survival amplitude for the state \& is given by 

S(t) = (^,e- lHt ^ 1 ^), (11) 

which is rigorously related to the energy distribution via 

~ iEt / n dW{E). 



S(t) = 

The survival probability 

P(t) = \S(t)f 



(12) 



(13) 



is thus given by Eq.(4). If, in a two body decay of a 
charged particle, dV(u>) is the probability of emitting soft 
photon radiation in the energy interval Tidu, then the 
radiatively corrected renormalization dWo(E) — » dW(E) 
is computed via the energy convolution 



( dW{E) x 
\ dE J 



'dWo(E-huj) 
v dE 



(14) 



If, during the decay, there are a mean number rik of pho- 
tons radiated into mode k with Poisson statistics, then 



dV(u) 
du 



£{n(- 



jS(u - y^rifc^fc). (15) 



{n} k k 

Employing the generating function in the time domain 



l dV{u) = e" x(t) , 



implies the simplification 

/"OO 

x(t) = (1 - 



(16) 



(17) 



where dN(u>) is the mean number of photons radiated 
into a bandwidth du>; i.e. Eqs.(15) and (16) imply 
Eq.(17) with 



dN(to) = ( nkS(uj — u>k) ) du>. 



(18) 



Furthermore, Eqs.(l),(3),(4),(14) and (16) imply 

P(t) = exp ( - Tt - 2 3?e x(t)) ■ (19) 

Eqs.(17) and (19) for the radiatively corrected survival 
probability are the central results of this section. Eqs.(7), 
(17) and (19) imply that the radiatively corrected transi- 
tion rate per unit time "/(t) as a function of time is related 
to the mean number of soft photons dN(uj) radiated into 
a bandwidth du via 



(20) 



j(t) = T + 2 / usin(ut)dN(u). 



The asymptotic Eq.(8) follows from Eqs.(9) and (20). 



III. COMPUTATION OF TRANSITION RATES 

In order to compute {dN {uj) / du) , and thereby 7(t), 
for the decay of a charged particle stopped in matter, we 
employ Schwinger's photon propagator method | pT| . For 
a charged particle moving on a path C, the self action 
due to virtual photons is given by 



2c 



T^nv{x - y)dx >1 dy v ', 



(21) 



c Jc 



where T> fJlV (x — y) is the photon propagator in the con- 
densed matter |12| wherein the original charged particle 
was stopped. Employing the A-space representation 



with 



L"(fc) = / e^dx", 



(22) 



(23) 



the action 



S = 



2c 



V 



(k)D^{h)L»{k)*(-^). (24) 



One may compute the mean number of radiated photons 
using 



N = 23m(S/ti). 
With the fine structure constant 
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(25) 



(26) 
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we have 



N = a 3m 



i[L^{k)D lxv {k)L v {k)*')( K 



( d A k 



(2tt)< 



(27) 



If Vi represents the initial four velocity of a particle (be- 
fore a two body decay) and Vf represents the final recoil 
four velocity of the produced charged particle after de- 
cay, then one easily obtains the usual expression [0 for 
L(k); i.e. 



^-'{(ri;) -(£;)}■ 

For the vacuum radiation distribution case 



r\Vac 



ac (k) _ (J^Jh^_\ 



(28) 



(29) 



so that Eqs.(37), (38) and (39) imply the well known g] 
vacuum radiated photon distribution 

*»-«-(^){(s^) -(&)}'■ <»' 

In the vacuum rest frame of the charged particle (before 
decay), we find the usual result, as in Eq.(18), 



fdio 



dN(Lu) = J 8(u- c|k|)d a iV, ac (k) = 0\^J . (31) 

For example, if in the rest frame of the charged particle 
(before decay), a single final state charge has velocity v, 
then 



«'>-(;){(£ WSS)-'}- 



Eqs.(20), (31) and (32) yield 
'2/?(v) 



j(t) = T + ^ — ^~~) ' ( vacuum decay). 



(32) 



(33) 



The radiative corrections for a charged particle decay, 
when the particle has been stopped in condensed matter, 
are somewhat more subtle. 



IV. DECAY OF PARTICLES STOPPED IN 
CONDENSED MATTER 

Condensed matter systems are often described (for 
complex frequency £ in the upper half 8m(() > plane) 
by a dielectric response function e(C), or by a conductiv- 
ity response function cr(C). These are related by 



<0 = i + 



(34) 



The above description may be incorporated into the tem- 
poral gauge photon propagator. £ — > (|u;| +i0 + ), written 
©as 

D(k,w) = 



-in 



|k| 2 - (u/c) 2 e(M 



i 1 C 2 e(|w|+iO+))}' 
(35) 

Employing the temporal gauge Eq.(35) in the evaluation 
of the number of radiated photons in Eqs.(27) and (28) 
yields 



(dN(bJ,v)\ 



(36) 



i.e. 



%,v)=(>/( 



'v ■ D(k,w + iO+)-v\ / d 3 k 



(k • v - uj) 2 



(37) 

The condensed matter version of the vacuum Eq.(32) is 
found (after some tedious integration) to be 

ne{ ( = )F(z(" + i0 + ), z> + *0+))} (38) 

where 

F(z,z*)=( 



(39) 
(40) 



and 



z + F 



(41) 



In the limit of real values for z 

'x + V 



YimT(x + iy) = {xln(^-—\ -2}, if \x\ > 1. (42) 
y^o I Vx — 1 / J 
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FIG. 1. Bremsstrahlung coupling strength as a function of 
velocity for an insulator from Eqs.(38) and (42). 
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In a (somewhat unphysical) model where the dielectric 
response s is independent of frequency, there will be a 
Bremsstrahlung regime for velocities obeying as shown 
in the above FIG.l; 



Thus, ha ~ 1 KeV which establishes the order of mag- 
nitude of the maximum frequency with which to define 
"soft photons" in this Ohm's law conducting model. 
Finally, for a conductor in the low frequency limit 



|v| < y~^/=J Brehmsstrahlung. 



(43) 



The value of the coupling strength in the insulating ma- 
terial Bremsstrahlung regime is given by 

«-."-G){(rW^H}- <«> 

In the high velocity regime there will be Cerenkov radi- 
ation. 



c > v > 



(£) 



Cerenkov. 



(45) 



However, the Cerenkov radiation regime can be discussed 
carefully only in models where the full complex dielec- 
tric response is taken into account. The dissipation in 
physical continuous media implies a finite /3(w,v) in all 
regimes. 

For example, in a model for an Ohm's law conducting 
material, the dielectric response function obeys 



, . , (Ama\ 



(46) 



For a non-rclativistic particle the resulting coupling 
strength /3(w,v) is plotted in FIG. 2. 



(£) - <' 



and 



8^ 



, as w^O. 



(47) 



(48) 



. From Eqs.(38), (40), (41), (47) and (48) it follows that 



A)(v)= lim/3(c,v) = 



, (conductor), (49) 



independent of the conductivity a. 
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FIG. 2. Coupling strength for a non-relativistic particle in 
a conductor. 

One may look at high velocity (for some values of 
(ui/a)) to note the almost discontinuous jump from the 
Bremsstrahlung regime to the Cerenkov regime in an 
Ohm's law metal. This is shown for (uj/cr) = 0.4 in 
FIG. 3. For a high conductivity metal, a <~ 10 18 /sec. 



FIG. 3. The almost discontinuous transition from the 
Bremsstrahlung to the Cerenkov regimes in a conductor. 



V. HOT SPOTS AND LONG TIME TAILS 

A typical model for the mean number of photons 
dN((jj) in a bandwidth dcu employs an exponential cut-off 
for high frequency; it reads 

dN(u;)=f3 a e-^(^-y (50) 

where (1/r) is the frequency cut-off. From Eqs.(20) and 
(50) it follows that 

7 (t) = r + 2/3 (^- I ) (51) 

may be used to compute the time dependent transition 
rate in Eq.(7). If t » t, then the time dependence of 
7(t) is given by Eq.(8). Ultimately lim t ^oo "f(t) = T, i.e. 



4 



the intrinsic decay rate. For very short times the intrinsic 
rate also dominates; i.e. linit—^o 7(i) = T. However, for 
intermediate times there exists a peak or "hot spot" in 
which the decay rate may increase substantially. The 
situation is shown in FIG. 4 which should hold true for 
both the vacuum decay and for decays in an insulator. 
For the insulating case j3o is plotted in FIG.l. 
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FIG. 4. Transition rate "hot spots" for various cut-offs. 



By plotting 7(t) in Eq.(7), one emphasizes the short 
time deviations from the uniform in time transition rate 
r. More conventionally H|, experimentalists di- 
rectly plot In P(t). In such plots, the short time hot 
spot appears not merely less pronounced, but in reality 
hardly visible. Thus, it may have escaped some deserved 
attention. 
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FIG. 5. In a direct plot of In P(t), the hot spot is not very 
prominent. 



From Eqs,(7) and (51), it is found in simple 
Bremsstrahlung models that 



P(t) 



t 2 



00 



exp 



(-rt) 



(52) 



To see what is involved, we plot the survival probability 
P(t) in the above FIG. 5. The short time hot spot for 
(Tt) 2 = 0.0001, which is so very obvious when plotted as 
in FIG. 4, is not at all so obvious when plotted as in Fig. 5. 
Both curves are mathematically equivalent in accordance 
with Eq.(7). To explore the short time hot spot on an ex- 
perimental level, one must examine the data in some de- 
tail during the period of the first life time of the survival 
probability, (say) at the early times < t < (0.5/r). For 
precision life time measurements, one tries a much wider 
time interval, (say) < t < (10/r). Such wide time in- 
tervals may mask important material effects which may 
have an effect on experimental precision. According to 
Eq.(52), the exponential decay law exp{— Tt) has a mate- 



when explored 



rials dependent prefactor (1 
over some twenty life times. 

Eq.(52) holds true for both conductors and insu- 
lators when the produced charged particle is in the 
Bremsstrahlung regime. The Cerenkov regime for insu- 
lators is a bit more subtle. If one may define for an 
insulator, a Debye relaxation time td, via 



3m(f(w + i0 + )) — > e((JTn), as uj — > 0, 



(53) 



then in the Cerenkov regime v > (c/^/e), and in the low 
frequency limit u> — ► 0, 

-(£)}■ (54) 

From Eqs.(20) and (54) it follows that in the insulating 
Cerenkov regime, the observed transition rate is renor- 
malized to 



(55) 



Finally, the case of an insulator with a fractal low fre- 
quency behavior flrjfl having exponent 77, 

5m(e(^ + z0 + )) -> g ( r~^ )» as w^0, (56) 



leads to 

7c(t,r?) =r + < <5 (7 ? )(^){ 



1 - 



sv 



OK 



t 



(57) 



where 



^(r 7 )=2r(7 ? )sin(^) ! rfo) = J°° x"e-*(fj). (58) 

From Eqs.(7) and (57) follows the Cerenkov fractal ex- 
ponent survival probability 
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P c (t)=exp(-rt-*c(*,»7,rD)) ) (59) 

where 

(60) 

and which exhibits (for short times) a stretched exponen- 
tial form. Note as r\ — > 0, Eq.(60) becomes equivalent to 
Eq.(55). 

VI. CONCLUDING REMARKS 

Deviations from the exponential laws in real time trig- 
gered by soft photons have been studied for decaying 
charged particles in materials. Under the standard as- 
sumption of factorizability of the "dynamical" energy 
distribution from the soft photon emission spectrum as 
given in Eq.(14), we have derived the transition rates for 
decays in insulating and conducting materials, both in 
the Bremsstrahlung and Cerenkov regimes. Several in- 
teresting results emerge some of which are listed below. 

If a particle decays in a conducting material, a novel 
discontinuity is predicted to occur. We show in Fig. 3, 
that as the velocity of the charged particle produced in 
the decay increases from low values (bremsstrahlung re- 
gion), the radiative coupling strength (and hence the 
transition rate) first decreases rapidly, shows a sharp 
discontinuity as it enters the Cerenkov region and then 
continues to decrease more slowly. Much care would be 
needed to experimentally observe such a discontinuity 
since it is rather sharp and hence confined to a very nar- 
row velocity range of the produced particle. 

Another prediction concerns "hot spots" and long time 
tails occurring both for the vacuum as well as an insula- 
tor. As exhibited in Fig. 4, the transition rate has a well 
defined maximum (hot spot) for intermediate times and 
is substantially different from its asymptotic value. On 
the other hand, the same effect is shown in Fig. 5 to be 
completely washed out in a standard logarithmic plot of 
the survival probability data commonly presented by ex- 
perimentalists. Thus, evidence, if any, for such an effect 
must be sought out through a careful study of the transi- 
tion rate in the short time interval, say < t < (0.5/r). 

The Cerenkov regime for insulators is also of particular 
interest since it leads to decays of the "stretched expo- 
nential" form. The fractal exponent in the absorptive 
part of the low frequency limit of the dielectric constant 
is shown to be directly related to the radiative exponent 
in the real time decay [see Eqs. (56-60)]. It would be 
worthwhile to check it experimentally. 

In view of the above predicted deviations from purely 
exponential decays, we urge that a concentrated, system- 
atic and precise experimental study of the transition rate 
be undertaken, both in the vacuum as well as in diverse 



materials for different decay particle speeds. Exponen- 
tial decays and Poisson statistics are almost axiomatic in 
experimental particle physics. Thus, any deviations are 
surely of fundamental interest. 
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